Flux-tunable barrier in proximity Josephson junctions 
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We report experiments on micron-scale normal metal loop connected by superconducting wires, 
where the sample geometry enables full modulation of the thermal activation barrier with applied 
magnetic flux, resembling a symmetric quantum interference device. We find that except a constant 
factor of five, the modulation of the barrier can be well fitted by the Ambegaokar-Halperin model 
for a resistively shunted junction, extended here to a proximity junction with flux-tunable coupling 
energy estimated using quasiclassical theory. This observation sheds light on the understanding of 
effect of thermal fluctuation in proximity junctions, while may also lead to an unprecedented level 
of control in quantum interference devices. 

PACS numbers: 74.45.-l-c, 74.50.-|-r, 74.40.-|-k, 74.78.Na 



The dynamics of a particle trapped in a shallow po- 
tential well is a fundamental problem that has wide ap- 
plicability to a number of areas in statistical physics. At 
finite temperatures, due to coupling to a thermal bath, 
the particle may escape from the potential well through 
the process of thermal activation over the barrier repre- 
sented by the edge of the potential. At low temperatures, 
there may not be enough thermal energy for the particle 
to overcome the barrier, but the particle may retain a 
finite probability of escaping from the well by tunneling 
through the barrier. 

For superconductors, two specific phenomena have 
been explored extensively in terms of the physics of ther- 
mal activation. The first phenomenon is the generation 
of phase-slips in a thin superconducting wire with cross- 
section size comparable to the superconducting coherence 
length ^5 P, In this case, thermally activated phase 
slips (TAPS) lead to the appearance of a resistance tail 
of the thin wire at temperature below the nominal transi- 
tion temperature Tc. The voltage V generated by TAPS 
is related to the time evolution of the macroscopic su- 
perconducting phase through the Josephson relation 
2eV/h = dip/dt, which on average is determined by num- 
ber of TAPS per unit time, each TAPS corresponding to 
a change of 27r. At the instant in space and time where 
such a phase slip event occurs, the superconducting order 
parameter vanishes, which costs an energy AF, the bar- 
rier over which the system must be thermally activated. 

The second phenomenon is the onset of finite voltage in 
a Josephson junction A resistively shunted junction 
(RSJ) can be modeled as a particle in a one-dimensional 
washboard potential in a viscous medium, where the dis- 
tance coordinate corresponds to the phase difference ip 
across the junction With no current through the 
junction, the system sits in a local minimum of the po- 
tential. Application of a current through the junction 
corresponds to tilting the washboard, shifting the posi- 
tion of the local minima. At some value of current less 
than the nominal critical current Ic, the washboard po- 



tential is tilted sufficiently for the system to be thermally 
activated over the barrier between two adjacent potential 
minima. Once this occurs, the system will continue to roll 
down the washboard potential, corresponding to a con- 
tinuous time evolution of ip, and hence a finite voltage 
will appear across the junction according to the Joseph- 
son relation. 

The Josephson junction case leads to an interesting ex- 
tension, where two junctions can be connected in parallel 
to form a dc superconducting quantum interference de- 
vice (SQUID). In this case, there are two independent pa- 
rameters (the phase differences across the two junctions), 
giving rise to a two-dimensional potential for the system. 
The system can transition from one local minimum to an- 
other through saddle-points in the potential 5]. Due to 
the fact that the position of the system on the two dimen- 
sional potential is sensitive to the external magnetic flux 
$, dc SQUIDs have been investigated extensively due to 
their device potential [1, 0|- With any real dc SQUID, 
the two Josephson junctions cannot be fabricated to be 
exactly the same, which restricts the ability to tune the 
system. Here we discuss doubly connected devices that 
are formed from normal metals in contact with supercon- 
ductors. The devices exhibit many of the properties of 
dc SQUIDs, but with the advantage that the devices can 
be designed to be almost perfectly symmetric, and hence 
allow unprecedented tunability of the system by means 
of an external magnetic flux. 

When a junction is made by a normal metal wire be- 
tween two superconductors, such an SNS junction ex- 
hibits properties different from a conventional tunnel- 
ing junction. The proximity to the superconductor has 
two major effects on the quasiparticles in the normal 
metal [8|. First, it induces superconducting- like corre- 
lations between quasiparticles that increase the conduc- 
tance of the normal metal, and second, it induces a gap 
in the quasiparticle density of states N{E) While 
in the superconductor the energy scale for the density of 
states is given by A, the energy scale for N{E) in the nor- 
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FIG. 1: (Color online) Schematic diagrams of two SNS quan- 
tum interference devices: asymmetric (a) and symmetric (b). 
The normal metal arms are shown in gold and the supercon- 
ducting wires are shown in gray. The arrow in both figures 
corresponds to the the direction of applied magnetic flux. In- 
set: Scanning electron micrographs of the devices measured, 
the scale bar is 1 fim. The calculated energy profile of the 
asymmetric/symmetric device are shown in (c)/(d) as a func- 
tion of the external magnetic fiux $ and the phase difference 
across the two superconductors (j>, which is determined by the 
external current through the device. The energy profiles are 
calculated based on the quasiclassical theory, as described in 
the text. 



mal metal is set by the Thouless energy Erh — hD/L?, 
where D = vpiji is the quasiparticle diffusion in the 
normal metal, i being the elastic scattering length, and 
L is the length of the normal metal. Similarly, the max- 
imal supercurrent that can flow through such a proxim- 
ity junction is set by Etk-, not by A for a conventional 
Josephson tunnel junction, when Etk ^ A (the long 
junction limit) [l^ . 

As with conventional junctions, two SNS junctions may 
be combined in parallel to form a dc SQUID. Schematic 
of two types of SNS junctions are shown in Fig. 1. In 
the asymmetric device shown in Fig. 1(a), the normal 
metal arms (shown in gold) connect to the superconduct- 
ing wires (shown in gray) at different points, resulting in 
likely different NS interface transparencies for the two 
arms of the device. Similar to a conventional dc SQUID, 
differences between the dimensions of the normal metal 
in the two arms, and differences between the NS inter- 
faces in the two arms of the loop invariably lead to an 
asymmetric dc SQUID [HI, However, the long-range 
nature of the Josephson coupling in SNS devices enables 
one to design and fabricate interference devices that are 
not possible with conventional tunnel junctions. One 
such device is shown schematically in Fig. 1(b), where 
the interference device consists of a single normal metal 
loop between the two superconducting contacts, so that 
the NS interface transparencies for the two arms of the 



device are the same. In this device, the modulation of 
quantum interference by an external magnetic flux occurs 
within the junction itself, i.e., the superconducting phase 
winding happens along the loop inside the junction, simi- 
lar to that in a superconduct ing loop that shows classical 
Little-Parks oscillations l^, 14, 1^. As we shall see, be- 
cause the NS interfaces are the same for both arms of 
the loop, this device behaves like a perfectly symmetric 
dc SQUID, while it is actually a single junction with a 
flux-tunable barrier. From another point of view, it can 
also be described as a thin su perc onducting wire with a 
flux-tunable phase-slip center [iq . 

To understand the modulation of thermal activation 
barrier in these SNS junctions, we calculate the energy of 
the system as a function of the phase difference Lp between 
the two superconductors . The energy of the system is 
given by [17 1 
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where Is{(p',^) is the supercurrent through the system, 
which is a periodic function of the phase difference if' 
and the externally applied flux $. To calculate the su- 
percurrent, we use the extended circuit theory [l^ and 
numerically solve the Usadel equations for the sample ge- 
ometries shown in Figs. 1(a) and 1(b). We assume here 
that the interface between the normal metal (N) and the 
superconductors (S) is perfectly transparent, the gap A 
regains its bulk value in the superconductor within a very 
short distance of the NS interface, and that the distribu- 
tion of quasiparticles in the superconducting reservoirs 
is given by the equilibrium Fermi function f{E). Since 
the characteristic unit for the supercurrent for an SNS 
junction is Exh/eR [l3|, we use {h/2e^R)Ec as the char- 
acteristic unit of energy for Ej. For the parameters used 
in this simulation, the amplitude of the supercurrent is 
about 0.2 Erh/eR, so the modulation of Ej at fixed $ 
is about 0.4 {h/2e'^R)Ec (note that since Ej has an arbi- 
trary constant from the integration in Eq. (1), in Fig. 1 
we assume Ej = a,t (p — —2tt). 

The resulting energy profiles for the system are shown 
in Figs. 1(c) and 1(d) as a function of ip and the two 
parameters under external experimental control. If 
is fixed at integral values of the superconducting flux 
quantum $0 = h/2e, for both geometries, there is an 
energy barrier for evolution of the phase, as shown by 
the trajectories of the blue particles at $ = 0, so that 
at low temperatures, the phase </) is stationary, and no 
voltage is developed across the device. For half-integral 
values of the apphed flux (<I> = [n -\- l/2)<i>o, where n 
is an integer) there is a difference between the asym- 
metric and symmetric cases. For the asymmetric case, 
there is still an energy barrier, as shown by the trajec- 
tory of the red particle in Fig. 1(c), although its height 
is smaller than that at integral values of the applied flux 
$ = n$o- Consequently, although the resistance of the 
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device at $ = (n + l/2)<I>o will remain finite at first even 
when the resistance at $ = n<i>o vanishes, the resistance 
at half-integral flux quanta will eventually vanish if the 
temperature is low enough, in the absence of quantum 
tunneling. In contrast, for the perfectly symmetric case, 
there is no energy barrier at half-integral flux quanta (as 
shown by the trajectory of the red particle in Fig 1(d)), 
so that the device will have a finite resistance even at the 
lowest temperatures, since Lp can evolve in time without 
any energy cost. 

Such a simplified physical picture can well explain the 
magnetoresistance of the SNS junctions shown the in- 
sets of Fig. 1(a) and 1(b). The details of fabrication and 
measurement are similar to those reported elsewhere [l^ . 
Figures 2(a) and 2(b) show the resistance of the asym- 
metric and symmetric devices respectively as a function 
of applied magnetic flux at a number of temperatures. 
The data are taken in the limit of zero dc current, with 
only a very small ac current for the resistance measure- 
ment (about 10-20 nA). For both geometries, the resis- 
tance is finite at all values of $ at higher temperatures, 
and is periodic in the applied flux, with a fundamental 
period of $o- As the temperature is lowered, the resis- 
tance for both devices vanishes near $ = n^o- As the 
temperature is lowered still further, the resistance of the 
asymmetric device also vanishes at half-integral values 
$=(71+ l/2)<i>o of the applied flux. In contrast, the 
resistance of the symmetric device at $ = (n -|- l/2)<i>o 
remains finite down to the lowest temperatures. 

The remarkable feature of the magnetoresistance 
curves for the symmetric sample is that peaks in resis- 
tance around <i> = (n-|-l/2)$o narrow as the temperature 
is lowered, saturating below about 0.15 K, as shown for 
the peak around $ = $o/2 in Fig. 2(c). According to the 
physical picture in Fig. 1(d), there is a small but finite 
barrier for the system to overcome, except at exactly half- 
integral values of flux. This barrier decreases monotoni- 
cally as the system approaches $ = (ri-|-l/2)$o- Thus, as 
the temperature is lowered, the system needs to be closer 
to <i> = (n -I- l/2)$o so that the particle can jump over 
the barrier and exhibit a finite resistance. At low enough 
temperatures, the system may tunnel through the bar- 
rier, resulting in a temperature independent resistance 
at lower temperatures. 

To quantitatively model the magnetoresistance around 
half-integral values of the applied flux for the symmetric 
device, we use the thermal activation theory proposed by 
Ambegaokar and Halperin (AH) for a resistively shunted 
junction as extended to a dc-SQUID. The AH the- 
ory predicts that the normalized resistance in the zero 
current limit is given by 



i?AH = V'(7/2) 



(2) 



where /q is the modified Bessel function, and 7 — 
2Ej/kBT is the ratio between the barrier and the ther- 
mal fluctuation energy. In the standard RSJ model, the 
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FIG. 2: (Color online) Resistance as a function of applied 
magnetic flux for the asymmetric sample (a) and the symmet- 
ric sample (b) at 0.6 K (blue curves), 0.4 K (green curves), 
and 0.03 K (red curves). While the oscillations for the asym- 
metric device die out rapidly with decreasing temperature, 
the oscillations for the symmetric device survive to the low- 
est measurement temperature, (c) Magnetoresistance of the 
symmetric device around $ = $o/2 at four different temper- 
atures. The solid lines are fits to the AH theory as described 
in the text. Below 0.15 K, the magnetoresistance does not 
change with temperature, (d) Critical current at zero applied 
flux for the symmetric and asymmetric devices. The plus 
symbols show the critical current expected from the fits of 
(c), multiplied by a factor of 5, as described in the text. 



Josephson coupling energy is Ej = {h/2e)Ic, Ic being the 
critical current. In our SNS junction case, Ej is given 
by Eq. ([T|), and is a function of the external flux $. As- 
suming symmetrical long SNS junctions, we obtain the 
usual sinusoidal dependence of the current Is on and 
the Josephson coupling energy 



E, 



2e '■ 



(3) 



where /,.($) — /c(0)| cos(7r<I>/$o)| as for a symmetric dc 
SQUID [l9|. Using this dependence for Ici.^), we have 
7 — hlc{^) /ekBT, and we can then fit the magnetoresis- 
tance of the symmetric device near $ = $o/2 using the 
AH theory with Eq. (2). 

The solid lines in Fig. 2(c) show the resulting fits to 
the AH theory at four different temperatures, using only 
the measured peak resistance i?p at $ = $o/2 and /c(0) 
as fitting parameters. The resulting values of /c(0) are 
a factor of 5 smaller than the experimentally measured 
values of /c(0) over the entire temperature range. Fig- 
ure 2(d) shows a comparison of /c(0) obtained from the 
fits (multiplied by a factor of 5) compared to the ex- 
perimentally measured values of /c(0). We note that the 
measured critical current saturates at lower temperatures 
and is about 10 times smaller than the value predicted 



4 



at zero temperature limit 1^, 11 1 , most probably due to 
an imperfect SN interface and heating [20|. 

The strong correlation between the measured and the 
fitted values indicates that the AH theory is applicable 
to diffusive SNS structures near $ = *&o/2, but the effec- 
tive potential barrier for thermal activation appears to 
be smaller by about a factor of 5 compared to Eq. ([3]). 
This discrepancy could not be due to quantum fluctua- 
tion [3] as this factor is temperature independent. It 
also can not be due to heating as then the measured crit- 
ical current and the activation barrier should be smaller 
than that inferred from the AH fit at zero current limit. 
In our opinion, this discrepancy could be due to the as- 
sumption of a standard external shunt resistor as the 
thermal noise source in the AH model, while the SNS 
junction is self-shunted. Another possibility is due to 
the neglect of inverse proximity effect and formation of 
minigap in both the normal metal and superconducting 
electrodes which may change the relation between 
supercurrent and energy in Eq. ([3|). 



We note in previous measurements |2l| . the current 
voltage characteristics of SNS junctions have shown large 
deviation from the AH theory and RSJ model, and theo- 
retically it is not clear how thermal fluctuation effects the 
proximity junctions. Here, being able to tune the barrier 
at zero current limit (close to equilibrium) enables us to 
apply the thermal activation theory to SNS devices for 
the first time. This finding provides a first step for a co- 
herent understanding of SNS devices and facilitates their 



applications such as sensitive detectors [22 



In summary, we have shown that diffusive SNS struc- 
tures have the potential to provide an unprecedented 
level of control in quantum interference devices. In par- 
ticular, we have demonstrated that one can make an SNS 
junction with a fiux-tunable barrier, which behaves like 
a perfectly symmetric dc SQUID. The potential barrier 
of such a device can be tuned to near zero with applied 
fiux, something that is extremely difficult to achieve with 
conventional dc SQUIDs. 
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